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Abstract 

We consider a large class of nearest neighbor attractive stochastic inter- 
acting systems that includes the asymmetric simple exclusion, zero range, 
bricklayers' and the symmetric K-exclusion processes. We provide exact 
formulas that connect particle flux (or surface growth) fluctuations to the 
two-point function of the process and to the motion of the second class 
particle. Such connections have only been available for simple exclusion 
where they were of great use in particle current fluctuation investigations. 

Keywords: simple exclusion, zero range, bricklayers', current fluctuations, sec- 
ond class particle, space-time covariance, two-point function, diffusivity 

MSC: 60K35, 82C41 

1 Introduction 

Serious research has recently been conducted on particle flux fluctuations in 
the asymmetric simple exclusion process (ASEP). This process is one of the 
simplest interacting particle systems, where particles are located at sites of Z, 
and each of them attempts to jump a unit step to the right after an independent 
exponential time with rate < p < 1, and to the left with rate q = 1 — p. If the 
destination site is occupied by another particle at the time of the attempt then 
the jump is suppressed. The translation-invariant extremal stationary measures 
for the ASEP are the Bernoulli distributions with density < Q < 1. 

In the Eulerian scaling, when time and space are rescaled with the same 
parameter, the hydrodynamic limit of the process leads to the inviscid Burgers 
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equation, see e.g. Rezakhanlou [T5] or Spohn [T7]. Burgers equation possesses a 
characteristic speed V e for each g value of the local Bernoulli equilibrium. 

Ferrari and Fontes 8J considered the time-integrated particle flux (or the net 
number of particles) that cross by time t the path of an observer who moves with 
a constant speed V. They proved that this random quantity is asymptotically 
normal and its variance scales with the first power of time when the observer's 
speed V differs from the characteristic speed V e . Fluctuations in this case 
basically come from the initial Bernoulli distribution. 

However, the correct scaling for integrated flux variance becomes t 2 ^ when 
the observer moves with the characteristic speed V = V s . This is the case when 
the genuine dynamical fluctuations of the process become visible. Many inter- 
esting results have been discovered in this direction, we refer to Prahofer and 
Spohn [12], Ferrari and Spohn [9], Balazs, Cator and Seppalalnen [3], Quastel 
and Valko [14], Balazs and Seppalalnen [5]. 

While the above results were true in simple exclusion, many other stochastic 
interacting systems are expected to have similar scaling properties. An exact 
connection between the flux variance and second class particles, that is derived 
with the help of the two-point function, proved to be essential in all of the above 
papers. In the present note we build on the first part of Balazs [2] to derive 
exact formulas that provide the same type of connection in a wide class of near- 
est neighbor attractive stochastic interacting systems. In [2] only asymptotic 
versions of these formulas were derived and used. The class we consider follows 
the basic ideas of Cocozza-Thivent [7], and it includes the asymmetric simple 
exclusion, zero range, bricklayers' processes, and the symmetric K-exclusion 
processes. Investigations on particle flux fluctuations will hopefully progress 
beyond ASEP to involve these more general models, and the present findings 
can serve as building blocks for such future arguments. 

Notation. Variables i, j, k and n will refer to sites of Z, while x, y and z 
will be integers in the single-site state space / to be defined below. 

1.1 The model 

The class of models described here is a generalization of the so-called misan- 
thrope process. We use a surface growth interpretation, but many members of 
this class can be understood in terms of particles jumping on the one dimen- 
sional lattice. For — oo < oj mm < and 1 < ui max < oo (possibly infinite valued) 
integers, we define the single-site state space 

I:={zeZ : LU min -1<Z< LU max + 1} 

and the configuration space 

n = {uj_ = {uoi) ieZ ■ Ui G 1} = I z . 

For each pair of neighboring sites i and i + 1 of Z, we can imagine a column 
built of bricks above the edge (i, i + 1). The height of this column is denoted 
by hi. A state configuration u £ 17 has components u>i — hi-i — hi El, being 
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the negative discrete gradients of the height of the "wall" . The evolution is 
described by jump processes. A brick can be added: 



or removed: 



with rate o^+i), 




(Wj, UJ i+1 ) ► (U>i + 1, LU i+l - 1) 

hi — > hi - 1 



with rate q(u>i, LUi 



+i) 



Conditionally on 0J.(t), these moves are independent. We impose the following 
assumptions on the rates: 

• The rates must satisfy 

(1.1) PK" 11 , • ) = P ( ■ , ^ laX ) EE ^ maX , • ) EE q( ■ , CJ™») EE 

whenever either ui min or aj max is finite. We assume that either p and q are 
non-zero in all other cases, or one of them is the identically zero function 
(totally asymmetric case). 

• The dynamics will have a smoothening effect when we assume monotonic- 
ity in the following way: 

, 12 > P(z + l,y)>p(z,y), p(y, z + l) <p(y, z) 

q(z + 1, y) < q(z, y), q(y, z + 1) > q(y, z) 

for y, z, z + 1 G I. This property has the natural interpretation that the 
higher neighbors a column has, the faster it grows and the slower it gets 
a brick removed. Our model is hence attractive. 

• We are going to use the product property of the model's translation- 
invariant stationary measure. For this reason, similarly to Cocozza-Thivent 
[7], we need two assumptions: 

— For any x, y, z £ / 

p(x, y)+p{y, z)+p(z, x) 
(1.3) +q(x, y) + q{y, z) + q(z, x)=p(x, z)+p(z, y)+p{y, x) 

+ q(x, z) + q(z, y) + q(y, x). 

— There are symmetric functions s p and s q , and a common function /, 
such that f(uj mm ) — whenever uj nun is finite, and for any y, z E I 
(1.4) 

p{y, z) = s p {y, z+ 1) ■ f(y) and q(y, z) = s q {y + 1, z) ■ f(z). 
Condition (|1.2p implies that / is non-decreasing on I. 
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• In order to properly construct the dynamics, restrictive growth conditions 
might be necessary on the rates p and q in case of an unbounded single- 
site state space /. We comment on this below. In particular, we assume 
that all moments of the growth rates are finite w.r.t. the distributions p,g 
introduced in Section [1.31 

At time t, the interface mentioned above is described by uJ_(t). Let ip : £1 — > K 
be a finite cylinder function i.e. ip depends on a finite number of ioi values. 
The growth of this interface is a Markov process, with the formal infinitesimal 
generator L: 



The construction of dynamics is available in the following situations. Several 
models with bounded rates are well understood and can be handled via the 
Hille-Yosida Theorem, see Liggett [11]. When the rates p and q grow at most 
linearly fast as functions of the local to values, then methods initiated by Liggett 
and Andjel lead to the construction of some zero range type systems (Andjel [JJ, 
Liggett [10], Booth and Quant [13]). The totally asymmetric zero range and 
bricklayers' processes with at most exponentially growing rates are constructed 
in Balazs, Rassoul-Agha, Seppalainen and Sethuraman [4]. See the definition of 
zero range and bricklayers' processes below. 

We assume that the existence of dynamics can be established on a set of 
tempered configurations (i.e. configurations obeying some restrictive growth 
conditions), and we have the usual properties of the semigroup and the generator 
acting on nice functions on this set. We also assume that f2 is of full measure 
w.r.t. the stationary measures defined in Section 11.31 Questions of existence of 
dynamics are not considered in the present paper. 

1.2 Examples 

There are three essentially different cases of these models. All of them are of 
nearest neighbor type. 

1. Generalized exclusion processes are described by our models in case 
both uj mm and w max are finite. 

• The asymmetric simple exclusion process introduced by F. 
Spitzer [TJ] is described this way by uj min = 0, w max = 1, f(z) = 
l{z = l}, 

s P (y, z ) =V Hv = z = 1} and s q (y, z) = q ■ l{y = z = 1}, 

where p and q are non- negative reals adding up to 1 (see (11.410 . In 
this case 

P(y, z ) =P'l{y = 1, z = 0} and q(y, z) = q-l{y = 0, z = 1}. 



U>j, ■ [ip{. . . , CJi - 1, UJ i+1 + 1, . . . ) - 



(1.5) 



+ 5> 



Wj+i) • [<p(. . . , u>i + 1, Wj+i - 1, • . . ) - <p(u)] . 
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Here u>i £ {0, 1} is the occupation number for site i, p{u>i, is 
the rate for a particle to jump from site i to i + 1, and q{u>i, Wj+i) is 
the rate for a particle to jump from site i + 1 to i. These rates have 
values p and q, respectively, whenever there is a particle to perform 
the above jumps, and there is no particle on the terminal site of the 
jumps. Conditions (|1.2p and (| 1 .3[) are also satisfied by these rates. 

• A particle-antiparticle exclusion process is also shown to de- 
monstrate the generality of the frame described above. Let p + q = 1, 

w mi„ = _ 1) ^max = L Fjx f(_^ = Qj f( Q j = Q ( creat i n), /(l) = O 

(annihilation) positive rates with c < a/2, 

s p (0, 1) =s p (l, 0) =p, s p (0, 0) = |^, 3 P (1, 1) = |, 
s,(0, 1) = 0) = g, s,(0, 0) = £, *,(1,1) = | 

and s p , s g zero in all other cases. These result in rates 

p(0, 0) = pc, p(0, -1) = p(l, 0) = ^, p(l, -1) = pa, 
9 (0,0)=gc, g(-l, 0) = g(Q, 1) = ^, g(-l, 1) = qa 

and zero in all other cases. If loi is the number of particles at site 
i, with u>i = — 1 meaning the presence of an antiparticle, then this 
model describes an asymmetric exclusion process of particles and 
antiparticles with annihilation and particle-antiparticle pair creation. 
These rates also satisfy our conditions. 

Other generalizations are possible allowing a bounded number of particles 
(or antiparticles) per site. 

2. Generalized misanthrope processes are obtained by choosing w mm > 



• The zero range process is included byp+g = 1, ui mia = 0, cj max = 
oo, an arbitrary nondecreasing function / : Z + — > R + such that 
/(0) = 0, 

Sp{y, z)=p and s q (y, z) = q, 

P(y, z) =p- f(y) and q(y, z) = q ■ f(z). 

Again, Ui represents the number of particles at site i. Depending 
on this number, a particle jumps from i to the right with rate p • 
f(ijJi), and to the left with rate g • /(wj). These rates trivially satisfy 
conditions ([L"2]) and (fO|l . 

3. General deposition processes are the type of these models where 
u> mm — —oo and o-i max = oo. In this case, the height difference between 
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columns next to each other can be arbitrary in Z. Hence the presence of 
antiparticles cannot be avoided when trying to give a particle representa- 
tion of the process. 

• Bricklayers' models. Let / : Z — > M + be non-decreasing, also 
having the property 

/0) • /(l - z) = 1 for all zeZ. 

The values of / for positive z's thus determine the values for non- 
positive z's. Define, with non-negative numbers p + q = 1, 

Spfai z) =p+ f and s g (y, z) = q + 



mm q mm 

which result in 

P{y, z ) =Pf(y)+Pf{-z) an d q{y, z) = qf{-y) + qf{z). 

This process can be represented by bricklayers standing at each site 
i, laying a brick on the column on their left with rate pf(—LOi) and 
laying a brick to their right with rate pf(u)i). They also remove 
a brick from their left with rate qf{uii) and from their right with 
rate qf(—u>i). This interpretation gives reason to call these models 
bricklayers' model. Conditions (| 1 . 2[) and (| 1 . 3[) hold for the rates. 



4. Symmetric processes are defined to have the identity p(y, z) = q(z, y). 
In this case (|1.3[) holds automatically, we only need to take care of (|1.2p 
and (fi~4l). 



• The symmetric K-exclusion process is obtained if we set a/ nm = 
0, uj max = K, f{z) = \{z > 0}, 

s p (y, z) = s q (y, z) = l{z, y < K}. 

These result in 

P(V, z) = q(z, y) = l{y > 0, z < K}. 

This process thus also has a product stationary distribution, as de- 
scribed below. 



1.3 Translation invariant stationary product distributions 

We now present some translation invariant stationary distributions for these 
processes. For many cases it has been proved that these are the only extremal 
translation-invariant stationary distributions. Following some ideas in Cocozza- 
Thivent [7j, we first consider the non-decreasing function / of (|1.4[) . For / 3 
z > we define 

z 

/(*)!:= n/0/)> 

y=l 
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while for I 3 z < let 

/W !: = V > 

n m 

y=z+l 

finally /(0)! : = 1. Then we have 

f(z)\ ■ f(z + 1) = f(z + 1)! 

for all z G /. Let 

5 f logfliminf (f(z)\) 1/z ) = lim log( f(z)) , if cj max = oo 

8 '. — \ V z — >oo / z — >oc 

oo , else 



and 



log ( limsup(/(-z)!) lA ) = lim log(/(-z)) , if u min = -oo 

— oo , else. 

By monotonicity of /, we have 9 > 9_. We assume 8 > 9_. With a generic real 
parameter 8 £ (d, 9) , we define the state sum 



e 

oo. 



Let the product-distribution ji have marginals 



1 e ez 

if z e /, 



^ e (z) = {w : Wi = z) : = { Z{9) f(z)\ 

if z £ I. 

Lemma 1.1. The product distribution is stationary for the process generated 
by (HU). 

We prove this lemma in Section [3] 



2 Results 

Introduce [a;] : = [x\ (floor) when x > 0, and [x] : = \x~\ (ceil) when x < 0. 
We start our process in the above described translation-invariant stationary 
distribution a The quantities of main interest are the following. For a fixed 

— a 

speed value V we define 

jW(t) :=h m (t)-h (0), 

the height of the column at site [Vt] at time t, relative to the initial height of 
the column at the origin. For V = we write 

(2.1) J(t) = J ( -°\t):=h (t)-h {0). 
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In particle notations of the models, J^ v \t) is the time-integrated particle flux, 
i.e. the net number of particles jumping through the moving window positioned 
at Vs, when s runs in the time interval [0, t}. Often we will make the choice 
MO) = 0. 

Theorem 2.1. For any V £ R and t > 0. 

oc 

(2.2) Var(J^(i))= ]T \[Vt] - n\ ■ Cov(w„(t), wb(0)). 

n=— co 

Moreover, we also have 

oo 

(2.3) n-Cov(o; n (t), w (0)) =*-Cov(p(o;o, wi) — g(wo, wi), (wo + wi)). 

n— — oo 

This theorem is proved in Sections [31 |3J and [5l Our primary objects of 
investigation are attractive systems. There is, however, only a minor point 
where our arguments use attractivity, and that is at the very end of Section 
16.21 where positivity of space-time correlations is used. This suggests that the 
above theorem should still hold in this class without attractivity, but we have 
not investigated this issue. However, the following results do genuinly require 
attractiveness. 

We introduce the notation 5^ £ fi, a configuration being one at site i and 
zero at all other sites. Let w £ f2 be such that ujq < cj max , and uj_ + : = lu + 5 Q . 
We say in this situation that we have a single second class particle between uj_ + 
and oj_ at site 0. Section HTT1 describes the basic coupling of two processes. With 
the above initial data this coupling conserves the single second class particle for 
all times t > 0: 

(2.4) u + (t)=u(t) + 6 Qit) . 

The quantity Q(t) is the position of the second class particle at time t, which 
performs a nearest neighbor walk, influenced by the ambient process w(-). 

We intend to consider the initial state uj in distribution (j, , together with 
a second class particle started from the origin. When u max < oo, this leads to 
a positive probability of ujq = cj max , in which case we cannot define our initial 
state lo + in il. We therefore introduce the indicator 

l{-}:=l{-}-l{^ <^ max }. 

With this notation it makes sense to have the 1-indicator of a second class 
particle-related event under an expectation. 

Theorem 2.2. We use the convention that the empty sum equals zero. Let 

g(z) := z ~^2y My)- 
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For any h£Z and t > 0, we have 




Moreover, the sum in the last display is non-negative for any loq . 

A short calculation shows that the mean of the sum on the right hand-side 
is Var(wo) (the variance w.r.t the distribution /ig). Denote by fl 9 the product 
measure of marginals fi g for all sites, except for the origin where 

(2.6) ]l g {uj_ : u = y}: = y , ; 2J 9(z)^e(z) (y E I). 

We write P and, correspondingly, E, Cov for probabilities of a process that is 
started in distribution 'p. . With this notation Theorem 12.21 rewrites as 

CovK(t), wo(0)) = Var(c^o) ■ P{Q(t) = n). 

Corollary 2.3. For any n S Z and ( e 1, i/ie variables wo(0) and w n (t) are 
non- negatively correlated. 

Corollary 2.4. Equation (|2.2I) can fee rewritten as 

(2.7) Var(j( y )(t)) = Var(w ) • E(|Q(t) - [Vt]|). 

Define now the quantity g = gg = Eg(u>o) which is a function of the pa- 
rameter 0. In particle systems this is simply the density of particles. Due to 
the definition of fig, we have dg/ d6 = Var(wo) > 0, which shows a one-to-one 
correspondence between g and 9. Define also the hydrodynamic flux TL{g) = 
Ee[p(wo, iO\) — q(u>o, u>i)] with the value of for which g = E(wo) = Eg^o) 
holds. The Eulerian scaling of these processes in many cases leads to Burgers- 
type hydrodynamic equations of the form 

d t g(T,X)+d x H(g(T, X)) = 0, 

where T and X are the rescaled time and space parameters (see e.g. Reza- 
khanlou [15] or Spohn [17]). It is well known that the characteristic speed of 
this PDE is V e = dH(g)/ dg which, by our definitions, can be shown to equal 
Cov(p(w , uji) - q(oj , Wi), (oJq + wi))/Var(w )- 

Corollary 2.5. With our new quantities, Equation (|2.3p takes the form 

(2.8) %{Q(t))=t-V e °. 

Remark 2.6. For the simple exclusion process u becomes the Bernoulli distri- 

V 

bution with density g — e e /(l + e e ), while J2 g is the same distribution conditioned 
on ujq = (which is the only way of initially having a second class particle at 
the origin). 
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Equation (|2.7p has appeared for the totally asymmetric simple exclusion 
process, with the special value V — V e (the characteristic speed), as Equation 
(1.7) in Ferrari and Fontes [S]. 

Corollary 12 . 51 has appeared in Prahofer and Spohn [TJ] for the totally asym- 
metric simple exclusion process. It is, on one hand, in accordance with the 
general phenomenon that the second class particle follows the characteristics of 
the hydrodynamic equation. On the other hand, J2 g is known not to be sta- 
tionary as seen by the second class particle in many models, which makes the 
validity of this finite-time formula somewhat surprising. 

Formulas (|2.7|) and (|2.8p have a great potential to serve as a basic building 
block in computing the order of flux fluctuations and diffusivity in our class of 
systems (see e.g. Quastel and Valko [14] for the definition of diffusivity in simple 
exclusion), as happened in the case of simple exclusion (Quastel and Valko [14], 
Balazs and Seppalainen [5]). 

3 Equilibrium and the reversed chain 

This section contains some basic computations regarding the measure p, . For 

w min _ 1 < Zj y < ^max + ^ define 

(3.1) p*(y, z) : = p(z,y), q*(y, z) := q(z, y) 
and 

(L*ij})(w) = ^p*(w l5 uj i+1 ) ■ [ip(. ..,Wi + l, u) i+ i - 1, . . . ) - ip(w)] 

(3.2) I£Z 

+ ^2q*(^i, • [ip(. ■ ■ , Wi - 1, u>i+i + 1, ...) - ^M] 

on bounded cylinder functions ip. 
Proposition 3.1. We have 

E • Lcp(u)) = E (<p(w) ■ L*ij){u)) 
on bounded cylinder functions ip and ip. 

Proof. By definition we have the property that for all uj aun — 1 < z < uj max , 

fxe(z) f{z + lY 
This, together with p.4[) , implies that whenever uj mm — 1 < y < uj max and 

P (y + l,z-l). ^^^ = P (z,y) 
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holds, and whenever cj min < y < uj max + 1 and u min - 1 < z < w max , then 

q(y -1,2+1) 7-r j— = q(z, y) 

holds. With p.ip . changing variables leads to 

( 3 - 3 ) E E p(l/.«)-G(y-l,«+l)w(y)w(«) 

y=OJ min 2— OJ min 
u; max u; max — 1 

= E E P(v> z)-G(y-l, z + l)fx e {y)fj. e (z) 

y — LU Tnin -\-l z=(jJ mm 

"ST* i , 1 -n ye(y + l ) Vei z - l ) n , \ , \ , \ 

= y. y. p{y + i,z-i) r-\ — n G w z)/w)M*0 

y— w min z— o; min + l 
(j max -l oj max 

= E E pi z ^y)- G iy-, z )My)^e{z) 

y=o; min z — LJ min + l 

w max w raax 



E E ■ z )My)w{z) 



y—UJ 1111 ' 1 z—iu" 



for any function G which makes the sums convergent. In a similar fashion, we 
get 



UJ UJ 



( 3 - 4 ) E E -z)-g(j/ + i, 



= E E li z , y) ■ G(y 7 z)n e (y)fig(z). 

y—u mlTL z — u; mln 

Let ^, (/? be bounded cylinder functions, and let I C Z be a finite discrete 
interval of which the size can be divided by three, and which contains the set 

{i G Z : if) ox if depends on k>j or on u>i-\} . 

Then the summation index i in the definition (|1.5|) of the generator can be run 
on the set X. We begin by changing variables u>i, Ui+i as in (|3.3|> and f)3.4j) : 



E (ip(cj) ■ L(p(w)) = E^p(w i+ i, Wi) • - •) ^ + 1= ^i+i - 1, •• • )vG*j) 

iei 

+ E}^q(oj i+1 , oji) ■ . ., 0Ji -1, +1, ... Mw) 

iez 
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Since |I| is divisible by three, (|1.3[) implies 

iex iei 
We thus conclude 

E(^(w)-i^) 

= E^p(wi + i, Wi) • [?/>(. .., Wi + 1, - 1, ...) - ^(w)]<p(w) 

+ E)]g(M( + i, Wi) • [ip(. . . , cot - 1, + 1, . . . ) - 

Comparing this display with (|3.2[) finishes the proof. □ 

Proof of Lemma \l.l[ The previous proposition with ^(w) = 1 shows that the 
expectation of the generator on any bounded cylinder <p is zero. □ 

Corollary 3.2. Formula ()3 . 2|) is the generator of the reversed process which 
has rates (|3.1[) . 

Note that the rates of the reversed process do not depend on the parameter 
9 of the original process' equilibrium distribution. 
Define the microscopic fluxes 

r(y, z) : = p(y, z) - q(y, z), r*(y, z) : = p*(y, z) - q*(y, z). 

The sum of the rates for a column will be 

S(y, z) : = p(y, z) + q(y, z). 

Corollary 3.3. 

E(r*(w , uJi) ■ (u>q - ui)) = -E(S(w , ui)). 

Proof. The proposition and Formula (|3.3p implies 

E(p* (w , wi) ■ (uj - ui)) = E(p(wi, w ) ■ {ujq - o>i)) 

= E(p(w , wi) • (wo - wi - 2)) 

= - E (P*( w o, ^i) • (w - wi)) - 2E(p(w , wi)), 

from which E(p*(wo, toi) ■ {ojo — u>i)) = — E(p(u)Q, 0Ji))- Similarly, 

~E(q*(uJo, lui) ■ (uq - ui)) = E(g(wi, uj ) ■ (w - Wi)) 

= E(g(w , w x ) ■ (w - wi + 2)) 

= -E(g*(w , wi) • (w - wi)) + 2E(9(w , Wi)), 

thus E(g*(w , wi) ■ (w - ^i)) = E(g(w , wi)). □ 
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4 Vertical growth 

It is easier to first consider the variance of J(t) (case V = 0). 
4.1 Martingale tricks 

For convenience, we introduce the notation A for the centered random variable 
A G L (/£.), and we further simplify notations by 

r(t) : = r(wo(*), wi(t)), r*(t) : = r*(w (t), wi(t)), S(t) : = S(w (t), w x (t)). 
Recall ([23]l . 
Lemma 4.1. 

t s 

Var { J (t)) =tE(S) + 2 J J E(r(«)r*(0)) cto ds. 
o o 

Proof. First notice that applying the generator on J k (t), then using Holder's 
inequality with the assumption that all moments of the rates p and q are finite 
implies in an inductive fashion that (the time-derivative of) any moment of J(t) 
is finite. By definition, E(J(i) | w(0)) = tr(Q) + o(t), hence 

i 

(4.1) M(t) :=J(t) - J r(s)ds 
is a martingale with M(0) = 0. Therefore, 

(4.2) Var( J(t)) = EM(<) 2 + 2 E^M(t) / r(s) ds) + e((/ r(s) ds 

o 

Due to E (M(t) 2 | w(0)) = 1 5(0) + o(t), the process 

t 

JV(t) : = M(t) 2 - [ 5(s) ds 



t 

2 







is also a martingale with -/V(O) = 0. Hence 

EM(t) 2 =tE(5). 

Using the martingale property of M, the second term of (|4.2|) can be written as 



t t 
2 y E(M(*)r(s)) ds = 2 f E(M(s)r(s)) ds. 
o o 
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Simply changing the limits of integration in the third term of (|4.2p , we have 

t t s 

E ((/ f ( s ) ds ) 2 ) = 2 / E ( F ( S ) J H u ) du ) ds- 



These calculations lead to 

t s 
(4.3) Var(J(t)) =tE(S*) + 2 J E(r(s)(Af(s) + / r(u)du))da 

o o 

t 

= tE(5) + 2y e(V(s) J(s)) ds. 
o 

In order to handle E (r(s) J(s)), we introduce *, the quantity corresponding 
to J in the reversed process by 

J (s) *(u) := J(s) - J(s-u) (s>u>0). 

This is the number of bricks removed from the column in the reversed process 
started from time s. As in case of J(t), a reversed martingale can be separated 

by 

u 

AfW*(u):= jW*(u)- [r*(s-v)dv. 



For this reversed object, AfW*(0) = and E (M< s )*(u) | = AfW*(s- 

£) if < s — £ < u, where stands for the natural filtration of the (forward) 
process. In view of this, 



E(r(s)J(s))=E[f( S )E(jW*( s )| 

= E I r(s) J r*{s-v)dv J = / E (r(u) r*(0)) dv, 

V o / o 

where in the last step we used time-invariance of the measure. Using this result, 
we obtain the statement from (|4.3[) by changing the order of integration. □ 

4.2 Space-time correlations 

In this subsection we denote r{u>i, and r(u>i, Wi+i) by r, and r,, respec- 

tively. We also keep the notation r — ro and r = Fo from the previous section. 
For jfe G Z, let 

d fc : f2 — ► 1 ; <4(u>) = ^fe 
be the fc-th coordinate of O. Then 

(4.4) (Ld fe ) (w) = r fe _i - r fe and (L"d*) (a;) = + r* k , 

where L* is the infinitesimal generator (|3.2p for the reversed process. 
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Lemma 4.2. For < a < 1 the functions 

oo oc 

(4.5) ip a :=^a k-1 d fc ^ a :=^a k d_ 



fc=i 



fc=0 



are /i^-a.s. well defined and finite, and 

\im(Lip a )(uj) = — \im (Lip a )(uj_) —r, 

a— >1 a— >1 

lim (L*^ Q )(a;) = — lim(L*<^ a )(a;)=r* 

a — >1 a — >1 

in £ 2 (^ e ). 

Proof. The a.s. existence of the sums above can easily be shown by using the 
Borel-Cantelli lemma for the sets 

A n := {lo : \u3 n \ > n}. 

We show the first equation for ip a . By (|4.4|) 

oo oo 

(4.6) (L<p a ) (w) = ro + (a - 1) nta*" 1 = r + (a - 1) ^ nta*" 1 . 



fe=X 



fe=i 



By independence of u>j and luj for i ^ j, E(r; • rj.) = if |Z — k\ > 1 and 



|E(rj • r fe )| = |Cov(n, r k )\ < Var(n) = E(r ; • n) 



12) 



if |fe — Z| = or 1. Hence the L 2 -norm of the second term on the right-hand side 
of (14.61) tends to zero as a — > 1: 



,fc-i 



fc=l 



<(a-lfJ2\\r k \\W k - 2 + 



k=l 



2(a-l) 2 ]T||n £ || 2 .a 



2fe-3 



fe=l 



(^-^11^(1 + 20^0. 
1 — ar a— >i 



The proof of the other three equations is similar. 



□ 



Lemma 4.3. Let $(0) = $(w(0)) 6e an L 2 -function that depends on the initial 
state only. Then 

t 

[ E(r(w) $(0)) dv - lim [E(£ a (i) $(0)) - E(£ a (0) $(0))] 



= lim [E(^ a (0) $(0)) - E(^ a (i) $(0))] . 
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Proof. We show the first equality, and by convenience we also center $(0) besides 



the centered rates. We first make use of the L convergence of the previous 
lemma: 

t t 

J E(r(v) $(0)) dv - limy E(L<p a (v) $(0)) dv 

o o 

t 

< limy \J E ( [r(v) - Ltp a (v)} 2 ) ■ E(5(0) 2 ) dv 



lim i-i/E (fr(Q) - Lip a (0)} 2 ) ■ E($(0) 2 ) = 0. 

a— >1 v 



Next we write 

t t 



E (5I a fc_1 £4(w(u)) 5(0)) dv. 
o fc=1 

Due to (|4.4p . a Cauchy inequality on E([rfc_i — r^] $(0)), finite moments of 
the rates and translation invariance, a k ^ 1 Ldk(u(v)) $(0) is absolute summable, 
"expectable" and integrable. Therefore 

t ^ t 

/oo _ 
E(L^ a (u) 5(0)) dw = Y,^ 1 / E(Ld fc (o;(t;))$(0))dt; 
£ — 1 

oo 

= ]T a*" 1 [E(d k (w(t)) $(0)) - E(4(#)) 5(0))] 
fc=i 

= E(<p a (t) 5(0)) -E(y> a (0)5(0)) 

by the integrated Kolmogorov equation on dk and another absolute integrability 
argument. □ 

Lemma 4.4. With $(0) as above, 
t 

E(F(-«) $(0)) d« = Hm [E$„(-t) $(0)) - E(VU0) $(0))] 
= lim [E(£ a (0) $(0)) - E(£ Q (-t) $(0))] . 

a — >1 

Proof. For the first line, repeat the previous proof except for the use of L* and 
tp a rather than L and ip a , and 

i 

y E(L*d_ fc (w(-v)) 5(0)) du = E(d_*(a(-*)) 5(0)) - E(d_ fe (o;(0)) 5(0)). 
o 

□ 
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Now we can compute the integrals in our expression for Var(J). 
Theorem 4.5. 

oc 

Var(J(i))= H-E(So(0)3„(t)). 



n— — oc 



Proof. The aim is to rewrite the double integral of Lemma 14.11 A slight mod- 
ification of Lemma 14.31 that includes a second integral as well implies that this 
can be done by 

Var(J(t)) =fE(S) + 2 lim / E(<p a (s) r*(Q)) ds - 2t Um E(£ a (0) r*(0)). 



The limits could be distributed to the difference since the second limit is finite 
by the centering and the product structure of /i^. The integral in this display 
is rewritten with a time-translation as 



t z 

J E($ a (s) r*(0)) ds = J E(£ a (0) ?*(-*)) ds 



= lim E(^ a (0) V 7 (-*)) - lim E(£„(0) V 7 (0)). 

7 — >1 7 — >1 

by Lemma 14.41 Notice again here that the second term on the right is finite. 
Hence with definitions (|4.5p , the variance of J(t) can now be written as 

Var(J(») = tE(S) - 2 Him E(ft,(0)r*(0)) 

+ 2 lim E(^ a (0)^~(-t))-2 lim E(^ a (0) ^L(O)) 

a, 7 — ► ! 0,7— »1 

= t E(5) - 2 1 lim E ^ a*- 1 Zo k (0) P (0) J 

(OO OC \ 

fc=i z=o / 

/ oo oo \ 

-2 lirr^E I ^a fc_1 a; A (0)^a'2_/(0) J . 



\k=l 1=0 



Using product property of the measure at time t = and the fact that r* 
depends only on luq and lo\ , most of our expressions become simple (recall that 
all quantities with tilde are centered random variables). We prove in Section 
that the limits in (14. T[) can be taken under the sum and the expectation: 



Var(J(i)) = tE(S) -2tE(2i(0)r*(0)) + 2^ j^E(2 fc (0)D_j(-t)) - 

k=l 1=0 

oo 

(4.8) =t~E{S) -2tE(r*(0)wi(0)) +2^nE(w n (t)u (0)). 



71=1 
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In the last step, we used space and time translation-invariance of the measure. 

We took advantage of the first equalities in both Lemma 14.31 and 14.41 The 
second identities therein can be used in a similar way to prove 

oo 

(4.9) Var(J(t)) = f E(S) + 2*E(r*(0) • S o (0)) + 2 ^ nE(S (0) £>_„(*)). 

n=l 

The statement now follows from Corollary 13.31 by taking the average of the 
previous two displays. □ 

Proof of (|2.3p . We only need to subtract (|4.9|) from (14. 8f> above to conclude 
the last statement of Theorem 12.11 Notice that r*(y, z) — r(z, y), but the 
symmetric term y + z makes the order of y and z in r(z, y) immaterial in the 
covariance. □ 



5 Non-vertical growth 

We turn to the variance of the quantity J^ v \t). First we consider V > values. 
Var(J^)(0)=E{(^ LV , J (t)) 2 } 

(5.1) = e{ ( v*j (*) - hvti (o)) 2 } - e{ (h [vti (o)) 2 } 

+ 2Ep LWJ (t)ft LvtJ (0)). 

Due to translation-invariance, the first term is Var(J(i)), computed in the pre- 
vious sections. By u>i — hi-\ — hi, ho(0) = and the product structure of the 
measure, the second term of the right-hand side of (|5.ip is 

(5.2) - E{ (h lvt} (0)) 2 } = - [Vt\ ■ E(£ (0) 2 ) = - [Vt\ ■ VarM. 

We compute the third term in the following lemma. 
Lemma 5.1. For V > 0, 

oo 

2E V*j(°)) = £ (\[Vt\-n\-\n\)-B(u n (t)oj o (0))+[Vt\-VaT(cJo). 

n— — oo 

Proof. Using tUi — hi-i — hi and ho(0) = again, 

yvt\ yvt\ yvt\ 

(5.3) E(V*j(*)v*j( Q )) = -E E ^ o (*)^( )) + EE E (^(*)^(°))- 

j=l i=l j=l 

With the martingale of (14. ip (for J(t) = h (tj), we have 

t 

E(/i (t)^(0)) = J-E(r Q (s)^(0))ds. 
o 
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We proceed via Lemma [ 

E(f (s)cj i (0))ds = lim E(p a (t)5,-(0)) - lim E(y> Q (0)G,-(0)) 



= ^E(^(*)^-(0))-VarK-(0)). 
»=i 

Again we justify in Section 16.21 that the limit can be taken under the sum on 
the space-time covariances. Plugging this in (|5.3[) gives 

oo yvt\ 

v(h vvt] {t)h\y n (<S)) =- H E (^(*)^(0))+ LVtJVar(wb) 

(5.4) I =L^tJ+U=i 

OO 

= -£](nA LV r iJ)E(« n (t)D (0)) + LVijVarfwo). 

n=l 

However, we also have, by the second identity of Lemma |4~31 
t 

/ E(r (s)£, (0))ds = - lim E(ip a (t) w,-(0)) + lim E(V> a (0) cj,-(0)) 





= - £ E(^(t)^(0)). 

i— — CJO 

Use this in (15.31) to obtain 



LvtJ vvt] 

^{h\yt\{t)h\y t ^))= J2 E E (^(*)^(°)) 
(5.5) l= -° OJ=1 



The statement now follows from taking the sum of (|5.4|) and (|5.5p . and from 
{\Vt\ -n+)-l{n < LVtJ}-(nA [y*J)-l{n > 1} = |L^J-n|-|n| (ra G Z). 

□ 



Proof of Theorem \2.1\ The case V" = and (|2.3|) are proved in the previous 
sections. For V > 0, combine (|5.1[) . Theorem 14.51 (|5.2[) . and the previous 
lemma. A similar computation shows the same result for negative V's. □ 



6 The second class particle 

In this section we show how to couple a pair of our models, with the help of the 
so-called second class particles. The space-time correlations seen in Theorem 
12.11 are rewritten in terms of the motion of this particle. 
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6.1 The basic coupling 



For configurations 77 and ( we say that 77 < (, if r/i < Q for all i 6 Z. Wc 
consider two realizations of a process, namely, ry(-) and £(■). We show the basic 
coupling which preserves 

(6.1) #)< C(*) 

if this property holds initially for t — 0. We say that di : = Q(t) — r/i(t) > is 
the number of second class particles present at site i at time t. 

The height of the column of £ (or rj) between sites i and i + 1 is denoted 
by 9i (or hi, respectively). Let gi f (or hi f) mean that the column of C (or 
the column of 77, respectively) between the sites i and i + 1 has grown by one 
brick. Similarly, gi J, (or hi j) means a brick-removal. Then the coupling rules 
are shown in Table [TJ Each line of this table represents a possible move, with 
rate written in the first column. As an illustration, we also indicate the change 
in the number di of second class particles at site i. These changes represent 
nearest neighbor walks of the second class particles, hence the total number of 
these particles is preserved. 



with rate 


9i 


hi 


di 


di+i 


p(Q, Ci+i) ~P(Vi, Ci+i) 


T 




I 


T 


pirn, V1+1) -p(m, d+i) 




T 


T 


I 


p(vi, d+i) 


T 


T 






q(0, 0+i) - Vi+i) 


i 




T 


i 


q(rn, rfi+i) - q(&, r) i+1 ) 




4 


I 


T 




I 


1 







Table 1: The basic coupling 

This coupling coincides with the well-known basic coupling for particle sys- 
tems. The rates of these steps are non-negative due to (|6.1I) and monotonicity 
(|1.2p of p and q. These rules clearly preserve property (|6.ip . since the rate of 
any move decreasing di becomes zero when di = 0. Summing marginally the 
rates of jumps of either rj(-) or £(•) shows that each process evolves according 
to its own rates. 

6.2 Space-time covariance and second class particles 

Recall the setting of Theorem 12.21 We now show how that theorem is derived 
from Theorem 12. II Recall that fJ,g(z) = for any z < uj mm or z > cj max . 

Lemma 6.1. For the pair (w(t), Q(t)) defined above Theorem \2.'2\ and for a 

function F : I — > R with F(uj max ) — and with finite expectation value 
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Y.F{z)^{z), 

"FMO)-l)/ieMO)-l) 



= E (1 



E(w n (t) : ^ UV " y .J'.^T'"' ---F"MO)) 

l{Q(t) = n}F(uJo(0))). 

Proof. Fix w min — 1 < z < aj max , and take conditional expectation of (|2.4|) : 

(6.2) E | wo(0) =z)=E (w„(t) | w (0) = z) + P (Q(t) = n | wo(0) = z) . 

Initially, w + (0) = uj(0) + S . Therefore, itself is also a process with initial 

distribution except for the origin. Hence 

E (co+(t) | lu o (0) = z) = E (w+(f) | a,+ (0) = z + l) = E (u„(t) | wo(0) = z + 1) , 

and (|6.2p can be written as 

E K(i) | w (0) = z + 1) - E (w n (t) | wo(0) = *) = P (Q(t) = n | uq(0) = z) . 

We multiply both sides with F(z) /ig(z) and then add up for all u; min — 1 < z < 
w™ to obtain 



J2 E K(t) | c^o(0) = z) • (F{z - 1) Me (z - 1) - F(z) /i (z)) 



(recall # e (z) = for z ^ 7 and F(w max ) = 0). The proof is finished by P(w (0) = 
z) = He(z). □ 



Proof of Theorem \2.2i By the previous lemma, our goal is now to find the cor- 
rect function F with finite mean, for which F(w max ) = and 



- F(z) = g(z) = z - ^2yn e (y) 



M*) yeI 
hold. By inverting the operation on the left side, we find 

Me (y) 



F(z)-.= ]T g(y) 



y= Z+ l WW' 

This function satisfies the conditions of the lemma, and (|2.5[) is proved. 

Now we turn to the proof that the sum in (|2.5[) is non-negative. Notice that 
E(g(ui)) = 0. The proof follows by both l{z > y} and g(z) being non-decreasing 
in z, and hence 

cj max 

< Cov(l{^ > y}, g(u)) = ^ ff(z)M*)- 

□ 
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We finally justify taking the limits under the summations in (|4.7p (and at 
some points in the proof of Lemma |5.1[) . For any a < 1 and 7 < 1, 



00 oc 00 oc 

^-^y E(|£ fc (0)| ■ ISM-*)!) <$> fc_1 E^ Var(o;o) < 00, 
fe=i ;=o fe=i ;=o 

and so we can reorder summations and expectations in (half of) (|4. T|) as 

00 oc 

£X> fc "V Cov(2 fc (0), 2_,(-t)). 
fe=l z=o 

Now the limits in a and 7 can be brought under the double sum by Corollary 
12.31 and Monotone Convergence. 
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